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We propose and experimentally realize a novel versatile protocol that allows the quantum state 
engineering of heralded optical coherent-state superpositions. This scheme relies on a two-mode 
squeezed state, linear mixing and a n-photon detection. It is optimally using expensive non-Gaussian 
resources to build up only the key non-Gaussian part of the targeted state. In the experimental 
case of a two-photon detection based on high-efficiency superconducting nanowire single-photon 
detectors, the freely propagating state exhibits a 67% fidelity with a squeezed even coherent-state 
superposition with a size |a|^=3. The demonstrated procedure and the achieved rate will facilitate 
the use of such superpositions in subsequent protocols, including fundamental tests and optical 
hybrid quantum information implementations. 

PACS numbers: 42.50.Dv, 03.65.Wj, 03.67.-a 

The optical hybrid approach to quantum information, 
which consists in combining the traditionally separated 
discrete- and continuous-variable states and operations, 
have seen important developments in recent years mm- 
First demonstrations of hybrid protocols offering novel 
capabilities have been achieved, including a long-distance 
single-photon entanglement witness based on local ho¬ 
modyne detections mm or the deterministic continuous- 
variable teleportation of discrete quantum bits [5] . Quan¬ 
tum repeater architectures have also been investigated 
laii] and the recent demonstrations of entanglement be¬ 
tween particle-like and wave-like qubits open the 

promise of heterogeneous networks based on the combi¬ 
nation of both encodings. This hybridization also finds a 
variety of unique applications in linear optical quantum 
computing mmm and quantum metrology m- 

Within this context, a considerable effort has been ded¬ 
icated to the generation of highly non-Gaussian states of 
light [2] . Specifically, free-propagating coherent-state su¬ 
perpositions (CSS), also referred as optical Schrodinger 
cat states, are an essential resource. Such states of the 
form |(a)±|—(a) consisting in a superposition of two coher¬ 
ent states with opposite phases and mean photon number 
\a/, plays the role of qubits in the coherent state basis 
[IMS]. The size \a/ of the superposition is a critical 
parameter as directly related to the overlap between the 
two coherent-state components. A value \a/ = 2 gives 
already an overlap \{a\ — c^)/ = 3.10“^. Be¬ 

yond their fundamental significance for hybrid Bell tests 
and the study of mascropicity for quantum optical states 
nanTj, the generation of CSS with this minimal size and 


a generation rate large enough to allow subsequent oper¬ 
ations will open a wealth of possible protocols and gate 
implementations [iiiaiiH]. However, such generation re¬ 
mains very challenging. 

In this endeavor, various optical circuits have been de¬ 
veloped to generate CSS using non-Gaussian resources. 
The first seminal scheme consisted in subtracting a 
single-photon from a single-mode squeezed vacuum ni- 
1^ . This operation can be performed by tapping a part 
of squeezed light and detect it with a single-photon detec¬ 
tor. This process results in heralding a squeezed single¬ 
photon, which exhibits a high fidelity with an odd CSS 
with \a/ rsu 1. In order to reach larger a values, other 
protocols have then been implemented. Two-photon sub¬ 
traction operated on squeezed light, with [22] or with¬ 
out [23] time separation, has led to values \a/ close to 
2. An alternative method has provided a 3 dB-squeezed 
CSS with \a/ = 2.6. For this purpose, a two-photon 
Fock state is first heralded, split on a beamsplitter and 
one output is measured by homodyne detection [24]. A 
quadrature measurement within a given acceptance win¬ 
dow is used as third conditioning. Recently, a similar 
approach based on the interference of two single-photons 
and homodyne conditioning has led to the same result 
[25] . In all these experiments, the generation rate is low, 
at the Hz level or below, precluding their use in sub¬ 
sequent protocols. This limited rate comes from either 
operation on single-mode squeezed light, where intrinsi¬ 
cally only a small part can be tapped, or the requirement 
of three cascaded conditioning procedures. 

In this Letter, we propose and implement a scheme en- 
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FIG. 1. (color online). Schemes for generating a|n —2) -\- I3\n) 
and theoretical simulations, (a) Two squeezed vacua are 
mixed on a tunable beamsplitter. The detection of n photons 
in one of the outputs heralds the generation. Equivalently, 
the scheme can be implemented directly from a two-mode 
squeezed vacuum. The two orthogonally-polarized modes are 
separated with a tunable mixing e induced by a rotated half¬ 
wave plate and a polarizing beamsplitter, with e = sin(2^). 
The angle ^ = 0° corresponds to the perfect separation, lead¬ 
ing to the generation of a n-photon state, (b) Two-photon 
detection: fidelity with a squeezed even CSS as a function of 
the ratio e/A. The fidelity is optimized by the size \a\‘^ and 
the squeezing of the target state. The inset gives the weights 
of the vacuum and two-photon components, (c) Three-photon 
detection and fidelity with a squeezed odd CSS. 


abling the heralded generation of coherent-state superpo¬ 
sitions with the minimal required resources. This state 
engineering relies on two-mode squeezed vacuum and a 
n-photon detection performed on one of the modes. No 
coherent displacement is required, in contrast to previ¬ 
ous works [261127]. The process allows to optimize the 
formation of the CSS in a versatile way, focusing all the 
non-Gaussian resources to prepare only the non-Gaussian 
part of the state. Experimentally, using a type-II optical 
parametric oscillator (OPO), we generate a state show¬ 
ing a 67% fidelity with a squeezed even CSS with a size 
\a\‘^=3. The preparation rate is two orders of magnitude 
larger than achieved heretofore with the aforementioned 
schemes. We also provide detailed characterization of the 
process enabling to change the size of the superposition 
and its squeezing. These observations are made possible 
by the combination of a large escape efficiency OPO and 
high-efficiency superconducting single-photon detectors. 

We first describe the proposed scheme, as illustrated on 
Fig. I^a). Two squeezed vacua, 5(^)|0), with a squeez¬ 
ing factor s = and A = tanh(^), are first overlapped 
with a 7r/2 phase-shift on a tunable beamsplitter. This 


configuration has been widely used in the balanced case 
as it provides a two-mode squeezed vacuum state for 
subsequent operations, such as continuous-variable tele¬ 
portation [28j |29]. In contrast, the beamsplitter here 
is asymmetric and a photon-counting measurement is 
performed on one of the outputs to herald the gener¬ 
ation of the synthesized superposition. In the follow¬ 
ing, we denote e ^ 1 the beamsplitter asymmetry, i.e. 
r = Y^(l + e)/2. In the ideal case of a photon-number 
resolving detector, the detection of n photons in the con¬ 
ditioning channel heralds the state 

2 ^f U\/nin-l)\n - 2) + A|n)) . 
yn(n — l)e^ + A^ ^ ^ 

This expression stays valid with multiplexed on-off single¬ 
photon detectors in the limit of low squeezing in order to 
neglect the multiple-photon components. 

Equivalently this scheme can be realized by 
starting from a two-mode squeezed state, i.e. 
photon-number correlated signal and idler beams 
\^)s,i = vTatvea”! n)s\n)i, as generated for 

instance by a type-II OPO below threshold. In this 
case, the modes are orthogonally polarized and the 
mixing is realized by separating the beams after a small 
polarization rotation induced by a half-wave plate. The 
mixing parameter is therefore given by e = sin(2^). The 
relative phase between the two CSS components can be 
easily controlled by birefringent elements inserted before 
the mode splitting. 

Numerical simulations are displayed in Figs. Sb) and 
(c) for a two- and a three-photon heralding detection. 
The fidelity is calculated with an ideal squeezed, odd or 
even, CSS as given by: 

5(^0|CSS±) = (|a) ± I - «)). 

For each point, the fidelity is optimized by the size and 
the squeezing factor s' = . Size \a\‘^ as large as 3 

for n = 2 and 5 for n = 3 can be obtained with fidelities 
above 98%. 

This procedure is an example of a very general method 
to optimize the preparation of non-Gaussian states with 
Gaussian operations [30]. This study proved that Gaus¬ 
sian operations, such as squeezing, can reduce the re¬ 
quired non-Gaussian resources for state preparation. By 
identifying the essential non-Gaussian operational cost, 
they defined a set of core non-Gaussian states that are 
good approximations of the targeted state via only a 
squeezing operation. This strategy can lead to better fi¬ 
delities with fewer photon detections than methods based 
only on state truncation. Within this approach, the work 
reported here provides the best approximation of such 
core states for the targeted CSS given available n-photon 
detection. In other words, it enables to produce the 
core state that requires the minimum number of pho¬ 
ton detections to be prepared. The remaining squeez- 
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ing operation can be implemented using for instance the 
universal squeezer already tested for non-Gaussian states 
m- Moreover, the core states are generally more robust 
against a damping than the superposition of coherent 
states itself [32]. They are therefore more suitable for 
the transmission and storage of non-Gaussian states. 

We now turn to the experimental realization. The 
setup is sketched on Fig. [^a). The two-mode squeezed 
vacuum is generated by a type-II OPO pumped far be¬ 
low threshold by a continuous-wave, frequency-doubled 
NdiYAG laser (Diabolo Innolight) [33]. The OPO is 
made of a triply-resonant semi-monolithic linear cavity 
where the input mirror is directly coated on one face of 
a 10-mm KTP crystal (Raicol) and the output coupler 
is a 38-mm-curved mirror. The input mirror exhibits an 
intensity reflection of 95% for the 532 nm pump and high 
reflection at 1064 nm, while the output coupler is highly 
reflective for the pump and has a 90% reflection for the 
infrared, leading to a bandwidth equal to 60 MHz and an 
escape efficiency estimated to ^ 0.9. 

At the OPO output, the orthogonally polarized signal 
and idler impinge on a half-wave plate and a polarizing 
beamsplitter, which enable to mix the two modes and 
separate the resulting orthogonal polarizations. An angle 
0 = 0° corresponds to separate signal and idler without 
coupling them. The reflected mode, which is used for 
the heralding, is frequency-filtered to remove the non¬ 
degenerate modes due to the OPO cavity. This filtering 
is based on a 0.5-nm interferential filter and a 320-MHz- 
bandwidth Fabry-Perot cavity [33l [34] . The filtered light 
is then split on a 50/50 fiber beamsplitter and detected 
by two superconducting nanowire single-photon detectors 
(SNSPDs) based on tungsten silicide (WSi) [35]. The 
optical stack of the SNSPDs was designed for maximum 
absorption at 1064 nm by optimizing the thickness of 
the dielectric layers. Details of the detectors system will 
be reported elsewhere. The system detection efficiency 
reaches 85% while dark count rate is below 10 cps, an 
important feature for high-fidelity state generation [36l 
[37]. 

As the experiment is performed in the continuous-wave 
regime, detection events can occur at different times. 
Here, the accepted coincidence window between the two 
heralding triggers is set to 0.8 ns, much smaller than the 
typical time given by the inverse of the OPO bandwidth. 
As a result, the temporal mode in which the state is gen¬ 
erated is given by a double-decaying exponential profile 
with a width at 1/e equal to 30 ns [38]. The heralded 
state is finally characterized by quantum state tomog¬ 
raphy via homodyne detection. Quadrature values from 
50000 measurements are processed with a maximum like¬ 
lihood algorithm [39], which provides the state density 
matrix and the corresponding Wigner function. 

The experiment is carried out in a cyclic fashion with 
50-ms data acquisition and 50-ms locking and calibration 
periods. During the locking time, two auxiliary beams. 


(a) 



FIG. 2 . (color online). Experimental setup, (a) A two-mode 
squeezed vacuum generated by a type-H OPO is split on a 
tunable polarizing beamsplitter. One of the outputs, which is 
frequency filtered by an interferential filter (IF) and a cavity 
(FP), impinges on a fiber beamsplitter. Coincidence events 
from high-efficiency superconducting single-photon detectors 
herald the preparation. The resulting state is characterized 
by homodyne detection, with an overall efficiency 77 = 85%. 
(b) Two-photon state generation when signal and idler are 
separated without mixing {0 = 0°). The first plot gives the 
measured quadrature distribution while the second one pro¬ 
vides the diagonal elements of the density matrix, with and 
without correction for detection losses. The right figure dis¬ 
plays the corresponding Wigner function without correction, 
with a negative ring value W = —0.13±0.01. The heralding 
rate reaches 200 Hz. 


otherwise blocked by mechanical shutters, are used. A 
first counter-propagating beam enables to lock the filter¬ 
ing cavity using a microcontroller-based technique m- 
A second weak beam is injected into the OPO cavity. 
This beam is amplified or de-amplified depending on the 
relative phase with the pump. To set this phase, the 
intensity is locked at a constant level. The interference 
between the seed and the local oscillator is then used 
to access the phase of the quadrature measured by the 
homodyne detection. 

As a first experimental characterization of the system, 
we consider the case 0 = 0°^ which should lead to the 
generation of a two-photon Fock state. Figure l^b) gives 
the measured quadrature distribution, the diagonal ele¬ 
ments of the reconstructed density matrix and the corre¬ 
sponding Wigner function. Without correction for losses, 
the two-photon component reaches a value of 58%. This 
value is the highest fidelity reported to date. As the OPO 
is pumped far below threshold, the three-photon compo- 
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FIG. 3. (color online). Experimental results. The plots give the density matrices (absolute values) and the associated Wigner 
functions for incremental values of the half-wave plate angle 0 (e between 0.017 and 0.16). The results are corrected for the 
r] = 85% detection efficiency. For 0 = 1.5°, the negative peaks of the Wigner function reach W = —0.27 ± 0.01. 


nent is kept below 3%. By correcting for detection losses, 
the two-photon component is 79 %, as mainly limited by 
the square of the OPO escape efficiency r]QpQ ^ 0.8. 

We now investigate the evolution of the heralded state 
when the signal and idler mixing is increased. Figure 
shows the density matrices and the corresponding Wigner 
functions for different values of the waveplate angle 
from 0.5° to 4.5°. As can be seen, the mixing strongly 
modifies the generated states, which are not anymore 
phase invariant and exhibit two negative peaks in the 
Wigner functions. These peaks correspond to oscillations 
in phase space between the two coherent-state compo¬ 
nents. Remarkably, due to the efficient SNSPDs and the 
strategy based on two-mode squeezed light, the herald¬ 
ing rate reaches 200 Hz. This rate is at least 2 orders of 
magnitude larger than the ones achieved heretofore in ex¬ 
periments usually based on photon subtraction operated 
on squeezed light or involving three conditioning steps. 

Turning to a characterization of the generated state, we 



FIG. 4. (color online). Fidelity between the generated state 
and a squeezed even GSS, for (a) ^ = 1.5° and (b) 0 = 2.5°. 
The plots give the calculated hdelity as a function of |a|^ and 
the squeezing in dB. The green crosses indicate the maximal 
hdelities. For 0 = 1.5°, the hdelity reaches 0.67 with an even 
GSS with a size \a\‘^ = 3 and a 4-dB squeezing. For 0 = 2.5°, 
the hdelity reaches 0.68 with an even GSS with a size |o|^ = 
1.9 and a 3-dB squeezing. 


present in Fig. [^its hdelity with an ideal squeezed even 
CSS, for two different angles. The hdelity is displayed as 
a function of the size \a\‘^ and its squeezing in dB. For 
0 = 1.5°, shown in Fig. [^a), a hdelity F = 0.67 ± 0.01 
is obtained for \a\‘^ = 3 and a squeezing of 4 dB. The 
generated state thereby exhibits the highest amplitude 
and hdelity reported to date for free-propagating CSS. 
Furthermore, the scheme is versatile. Indeed, by simply 
adjusting the mixing ratio e, i.e. the wave-plate angle, the 
properties of the state can be engineered. For instance, 
as seen in Fig. Sb), the state generated with 0 = 2.5° 
exhibits a maximal hdelity F = 0.68 ±0.01 for \a\‘^ = 1.9 
and a squeezing of 3 dB. For any targeted squeezing, one 
can hnd the optimal angle to be used. 

We hnally note that, given the escape efficiency of the 
OPO, the corrected hdelity is expected to reach a value 
close to 80%, similarly to the one obtained for the two- 
photon state displayed in Fig. ib). The discrepancy 
mainly comes from the strong phase sensitivity of such 
superposed states. Due to phase calibration performed 
only every few seconds at the beginning of each acqui¬ 
sition sequence and transient perturbation of the exper¬ 
iment by mechanical shutters, errors on the phase cali¬ 
bration reduce the achieved hdelity. 

In conclusion, we have proposed and implemented a 
scheme enabling a versatile quantum state engineering 
of squeezed coherent-state superpositions, which are 
essential highly non-classical states in the emerging 
held of optical hybrid quantum information processing. 
Given a n-photon detection, the protocol provides the 
best approximation of the non-Gaussian core state for 
the targeted CSS. Experimentally, in the case of a 
two-photon heralding, a 4 dB-squeezed even CSS state 
is obtained with \a\‘^ = 3 in a well-controlled spatio- 
temporal mode. Importantly, this scheme is based on 
twin beams, and not on single-mode squeezed light or 
cascaded conditionings that usually strongly limit the 
heralding rate. The achieved preparation rate makes 
these states immediately suitable for subsequent experi- 
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ments where they can be used as initial resources, such 
as in test implementations of all-optical non-Gaussian 
gates for quantum communication networks [HITIITH]. 
From this perspective, it is also advantageous that 
the generated squeezed CSS can be optimally adapted 
for transmission in optical channels and storage in 
atomic memories [32]. Further improvements of the 
non-Gaussian core states can be achieved by feasible 
extensions with developing photon-number resolving 
detectors. Beyond direct applications, the generated 
CSS and its extension to larger core states can also be 
used to investigate processes connected to fundamental 
quantum-classical transition mi- 
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